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1. Introduction . Present quantum theory does not make definite prediction of the value of an 
observable in an individual observation except in an eigenstate of the observable. Application of 
quantum rules to two separated systems which interacted in the past together with a local reality 
principle (Einstein locality) led Einstein, Podolsky and Rosen^ to conclude that quantum theory 
is incomplete. BelP showed that previous proofs of impossibility of a theory more complete 
than quantum mechanics^ made unreasonable assumptions; he went on however to prove^ that 
a hidden variable theory agreeing with the statistical predictions of quantum theory cannot 
obey Einstein locality. 

Bell's research was influenced by the construction by De Broglie and Bohm^ (dBB) of a 
hidden variable theory which reproduced the position probability density of quantum mechanics 
but violated Einstein locality for many particle systems. For a single particle moving in one 
dimension with Hamiltonian 

H = -h^/{2m)d^/dx' + U{x), (1) 

and wave function ilj{x,t), de Broglie-Bohm proposed the complete description of the state to 
be {A(t), where X{t) is the instantaneous position of the particle, and its momentum is 

PdBBiX, t) = md\/dt = [Re ri-ifi ^^l^/^x)/{\^l;\^)l=x■ (2) 

In an ensemble the position density p(A, t) agrees with \ i/j{X,t)\'^ for all time. However, Takabayasi^ 
pointed out that the joint probability distribution for position and momentum given by the 
theory 

PdBB(A,p,i) = \il{X,t)\'^6 {p - PdBB{\t)) (3) 

does not yield the correct quantum mechanical expectation value of for integral n 1. 
De Broglie^ stated that these values in his theory "correspond to the unobservablc probability 
distribution existing prior to any measurement" and measurement will reveal different values 
distributed according to standard statistical quantum mechanical formula. On the other hand 
position measurements have no central role since they simply reveal the existing position. The 
asymmetrical treatment of position and momentum in the dBB theory constitutes breaking of 
a fundamental symmetry of the quantum theory and has been considered by some physicists 
as a defect of the dBB theory (Holland, Ref. 5, p. 21). 
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Without using hidden variables, Griffiths'^ and Gell-Mann and Hartle* introduced joint 
probability distributions for noncommuting observables at different times in the consistent 
history approach to quantum theory of closed systems. Wigner^ had earlier introduced a joint 
distribution for x and p at the same time, 

/oo (j^y / y \ f y \ 
oo 27r^^* \ ^ 2 ' 7 ^ V ~ 2 ' 7 ^^P(^^'^/^) (^) 

which yielded the correct quantum probability distributions separately for x and p on inte- 
gration over p and x respectively. The Wigner distribution cannot however be considered a 
probability distribution because it is not positive definite, as seen from the fact that the integral 

j dx dp PwA^,p)pwA^^P) = KV^^IVIStt^) 
vanishes for two orthogonal states ip, 0. 

We wish now to propose a deterministic quantum theory of a closed system with the fol- 
lowing properties. (We consider in this paper only 1 particle in 1 space dimension). 

(i) The system point {x{t),p{t)) in phase space has a Hamiltonian flow with a c-number causal 
Hamiltonian Hc{x,p, ip{x, t),t) so that in an ensemble of mental copies of the system the phase 
space density p{x,p,t) obeys Liouville's theorem 

dp{x,p,t)/dt^O. (5) 

Here ip{x,t) is the solution of the usual Schrodinger equation 

ihdi;{x,t)/dt = H^{x,t) (6) 

with H being the standard quantum mechanical Hamiltonian for the system and He being 
determined from the following criteria. 

(ii) Each pure "causal state", i.e., a set of phase space points moving according to a single 
causal Hamiltonian He has phase space density of the deterministic form 

p{x,p,t)^\iP{x,t)\^5{p-p{x,t)), (7) 

in which p — p{x, t) — not only determines p as a function of x, but also determines x as a 
function of p at each time (step functions being allowed when necessary). Eqn. (7) guarantees 
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on integration over p the correct quantum probability distribution in x for any real function 
p{x,t). The function is determined from the requirement that on integration over x, p{x,p,t) 
should also yield the correct quantum probability distribution in p. That such a determination 
is possible and unique apart from a discrete 2-fold ambiguity will be a crucial part of the present 
theory. It is obvious that our p{x, t) will have to be different from the p{x, t) of de Broglie-Bohm 
theory. 

(iii) Since the quantum probability distributions for x and p in the statistics of many measure- 
ments are exactly reproduced, so are the standard uncertainty relations. 

In Sees. II, III we describe the construction of the momentum p{x, t) and the causal Hamil- 
tonian He, in Sec. IV applications to simple quantum systems, and in Sec. V conceptual 
features of the new mechanics. 

2. Construction of Joint Probability Distribution of position and 

momentum. We seek a positive definite distribution of the form (7) where p is a monotonic 
function of x 



The monotonicity property ensures that for a given t, the 5-function establishes one-to-one 
invertible correspondence between x and p whenever dp/dx is finite and non-zero. (This is 
the simplest qualitative assumption about p{x, t) which will be shown to result in Hamiltonian 
evolution; in future development we should try to replace the monotonicity assumption by the 
assumption of Hamiltonian evolution). The requirement of reproducing the correct quantum 
probability distribution of p is that 



e dp{x,t)/dx >0, e — ±l 



(8) 




(9) 



where ip{k,t) is the Fourier transform of ip{x,t). We substitute the ansatz (7) into (9) and 



integrate over momentum to obtain 



J —oo Jp\ 



p{x',t)<p 



dx'\i^{x',t)fS{p'-p{x',t)) = £^^|V5 (pt 



(10) 
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The region p{x', t) < p becomes x' < x if e = 1, and x' > x if e = —1, where p{x, t) — p. Thus, 
we obtain, for e — ±1, 



The left-hand side is a monotonia function of x which tends to 1 for ex — > oo for a normahzed 
wave function; the right-hand side is a monotonia function of p tending to 1 for p ^ oo 
(Parseval's theorem). Hence, for each t, Eq. (11) determines two monotonia functions p of x, 
one for each sign of e. (Note that the curve p{x, t) may have segments parallel to x-axis or p-axis 



so determined yield via Eq. (7) phase space densities p±, with different causal Hamiltonians 
{Hc)± determined below. 

3. Determination of the Causal Hamiltonian. We view p{x,p, t) as describing an ensem- 
ble of system trajectories in the phase space. We saw in the last section that such a description 
is possible at each time. We would now like to find causal Hamiltonian such that the time 
evolution in phase space implied thereby is consistent with the time dependent Schrodinger 
equation. 

In order that the total number of trajectories is conserved in time we must have the conti- 
nuity equation 




(11) 



corresponding to ■0(x, t) or ip{p/h, t) vanishing in some segment). The two curves p — p±{x, t) 



dp/dt + d{px) /dx + d{pp) /dp = 



(12) 



If the dynamics of the trajectories is of Hamiltonian nature i.e. 



X = dHc/dp, p 



dHc I dx 



(13) 



then we have Liouville's theorem that the phase space density is conserved. 



dp/ dt -\- xdp/ dx -\- pdp/ dp = 



(14) 



I.e. 



dp/dt + (dHc/dp) dp/dx — {dHc/dx) dp/ dp — 0. 



(15) 



4 



The c-number Hamiltonian He describing the causal time evolution of the trajectories in the 
phase space will be allowed to be different from the usual ^-number Hamiltonian H describing 
the time evolution of the Schrodinger wave function according to Eq. (6). 

On substituting into Eq. (15) the ansatz (7) discussed in the last section, we obtain 

iKp-p) + ^mp-p))-^ (16) 

where 

^ = d\ilj\^/dt + (dHc/dp) d\^\'^/dx - dr]/dp 

Tj = — IV'P {dp/dt + (dp/dx) dHc/dp + dHc/dx} . 

We thus need for consistency 

^ = and 77 = Oifp=p. (17) 

We now specialise to the usual case when H is given by Eq. (1). We find that this situation is 
taken care of with the choice of Hc{x,p, t) 

Hc^^{p-A{x,t)y + V{x,t). (18) 

The causal Hamiltonian is of the Newtonian form apart from the introduction of a vector 

potential A{x,t) and allowing the potential V{x,t) to differ from U{x). Eqs. (17) lead to 
the following equations to determine V and A (after using Schrodinger eqn. to substitute for 

d\^\ydt), 

-dV(x, t) /dx = dp(x, t) /dt + (2m)-^d (p(x, t) - A(x, t)f /dx, (19) 
d[\ij\'^{p- A-mv)]/dx^O, (20) 

where v is given by 

v{x,t) = h/{2im) den{ilj/^*)/dx (21) 

which is just the de Broglie-Bohm velocity. Eq. (20) implies that the quantity in square 
brackets must be a function of t alone. We choose this function of t to be zero in order to avoid 
a singularity of the vector potential at the nodes of the wave function. We thus obtain 

A{x,t) ^p{x,t) -mv{x,t) (22) 
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With the calculation of the causal Hamiltonian thus completed via Eqs. (18), (19) and (22) a 
consistent Liouville description emerges. 

It should be stressed that the qualitatively new feature of the theory, p{x,t) ^ mv{x,t) is 
quite independent of the specific ansatz (18) for He- Comparison of the continuity equation 
for the spatial probability density p(a;, t) in a deterministic theory with that following from 
Schrodinger eqn. plus the requirement p{x,t) = \ip{x,t)\'^ imply that dx/dt = v{x,t), the dBB 
velocity. The assumption p — mv will then lead to the dBB answer for momentum probability 
density which conflicts with the quantum answer. Hence, to reproduce both position and 
momentum probability densities correctly we need p ^ mv. On taking ensemble average, the 
equality is restored in our theory in agreement with Ehrenfest's theorem. 



4. Illustrative Examples, (i) Quantum Free Particle. Let the quantum free particle be 
described by the Gaussian momentum space wave function 



i^ip/h,t) = (Q;7r)-^/^exp f-(p - (3f/{2ah^) - ipH/{2mh) 



so that the coordinate space wave function is 



xp{x, t) — (ttq;) ^^'^ {ma/ (m + iaht))^^^ exp /, 



/ = - W2) 
Our procedure yields 



{x — (5t/mY — i 



p-/3 = ±h\ 



' aht 



-x^ + 



2/3x pH 



m 



m^a^ 



ah mah^ 



1 + 



aW 



m^ 



A={p-P) 



m^ + a'^hH'^ 

hat 



X 



m 



It 



,Jm'^+{hatf_ 



and 



dV/dx = ±(m^ + a'^hH^)-^''^[xt{ahf + (5m]{ham) 



(23) 



(24) 



(25) 



(26) 



The determination of the causal Hamiltonian is now complete apart from an irrelevant additive 
function of t. The quantum potentials A and V are seen to be proportional to h in this example. 
An interesting feature of Eq. (25) is that for e = +1, for t » m/ {ah) fa 2{hE/ {AEy), it agrees 
with the naive classical expectation corresponding to zero vector potential. 



(ii) Quantum Oscillator. For the minimum uncertainty coherent state of the harmonic oscillator 
of mss m, angular frequency cu and amplitude of oscillation a we find 



. . muj 
p{x,p,t) = ^/ — exp 



[x — acos[Oi)t)) 



S{p-p{x,t)), (27) 



where 



and 



p{x,t) — —mwasm{u!t) ± mw{x — acos{u!t)), (28) 



A{x,t) — ±mu;{x — acos{u;t)), (29) 
—dV{x, t)/dx = —muj'^a cos(a;t) ± miu^a sin(a;t) (30) 
The causal Hamiltonian yields the equation of motion 

md^x/dt^ — —muj^acosut (31) 

which results in exact harmonic motion even for x away from the centre of the packet. We do 
not of course expect this for solutions of the Schrodinger eqn. different from the coherent state 
here considered. 

5. New conceptual features, (a) We have derived corresponding to every quantum wave 
function V', two joint probability distributions for position and momentum of the form (7) which 
are (i) positive definite, (ii) have Hamiltonian evolution with causal Hamiltonians {Hc)± and 
obey (iii) 

/ {f{x)+g{p))p±{x,p,t)dx dp = (^j, (^f{x)+g(^-ih^'^ , (32) 

for arbitrary functions f{x) and g{p). Eq. (32) is the major advantage of the present theory over 
the dBB theory. We postpone the discussion of measurements until we present a generalization 
of the theory to many particles. 

(b) Since both p+ and p_ obey Eq. (32) so will p — Cp+ + (1 — C)p_ with < C < 1. But since 
p+ and p_ correspond to different causal Hamiltonians {Hc)±, p will not correspond to a 'pure 

causal state'. We are led to the concept of a pure causal state as being more fine grained than 
a pure wave function ip. All p — Cp+ + {1 — C)p- correspond to ip {p ip) for a continuum 
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of values of C, but only C = 0, 1 correspond to pure causal states. To quantum density matrix 
states TiCa\il^a){i^a\ correspond phase space densities YlCaPa if Pa ^ i^a- 

(c) One can ask if Eq. (32) can be generalized to more general quantum obscrvablcs. Here we 



face the old problem that there exist nonclassical observables e.g. x i—ih4-\ x, ((—ih^j xxh.c.) /2 



which have different expectation values but the same 'naive' classical analogue x^p. A trivial 
but nonunique way followed already for the dBB distribution is: given a nonclassical observable 
A the phase space analogue can be f{x,p,ip) such that f{x,p,ip) — ip'^Aip/\ip\'^. 

(d) It can be proved that for the individual trajectories, Newton's first law [dp/dt = for 
U (x) — 0] holds, unlike dBB theory. 

(e) A qualitative advantage of our theory over the dBB theory is the symmetric treatment of 
X and p obvious from our phase space density 
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higher dimensions. We thank C.V.K. Baba, R. Cowsik, D. Dhar, P. Eberhard, B. d'Espagnat, 
J. Finkelstein, P. Holland, D. Home, S.S. Jha, L. Kadanoff, A. Kumar, P. Mitra, G. Ra- 
jasekaran, N.F. Ramsey, D. Sahoo and E.J. Squires for discussions and correspondence on an 
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p{x,p,t) = mx,m^{p/h,t)\'h-'S[ JZ^dx'\^{ex',t) 
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